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Abstract 

O 

We discuss the possible role of quantum horizon fluctuations on black hole radiance, especially 

(N 

%4 . whether they can invalidate Hawking's analysis based upon transplanckian modes. We are partic- 

ularly concerned with "enhanced" fluctuations produced by gravitons or matter fields in squeezed 
vacuum states sent into the black hole after the collapse process. This allows for the possibility of 
increasing the fluctuations well above the vacuum level. We find that these enhanced fluctuations 

O'. 

could significantly alter stimulated emission, but have little effect upon the spontaneous emission. 
Thus the thermal character of the Hawking radiation is remarkably robust. 



> 

o 

OO 
On 

cn 

o 

OO 
O 



X 



PACS numbers: 04.70. Dy, 04.60.Bc, 04.62,+v 



*Electronic address: robert@cosmos.phy.tufts.edu 
^Electronic address: ford@cosmos.phy.tufts.edu 



1 



I. INTRODUCTION 



Hawking's discovery [l| of black hole radiance has forged an elegant link between rela- 
tivity, quantum theory, and thermodynamics. However, some unsolved problems remain, 
including the information and transplanckian issues. The question of whether information 
is lost in the black hole evaporation process has been vigorously debated by many authors. 
(See, for example Ref. [2(] and references therein.) The transplanckian issue arises because 
Hawking's original derivation used quantum field theory on a fixed background spacetime 
and requires incoming vacuum modes with frequencies far above the Planck scale. Alter- 
native derivations have been proposed, especially by Unruh [!, H[ and by Jacobson and 
coworkers [1, @, 0, @, E| which involve cutoffs and a non-linear dispersion relation. The non- 
linearity can lead the the phenomenon of "mode regeneration" , whereby the modes needed 
for the outgoing particles are created just before they are needed, rather than being red- 
shifted transplanckian modes. However, this approach requires a breakdown of local Lorentz 
invariance and hence new, as yet unobserved, physics. 

In this paper, we wish to consider the effects of quantum horizon fluctuations on the 
Hawking process. This is a topic which has been discussed by several authors from various 
viewpoints. An early discussion was given by York [l(|, who used a semiclassical approach 
to estimate the magnitude of the quantum metric fluctuations near the horizon. Ford and 
Svaiter [Hj used York's estimate to treat fluctuations of the outgoing rays, and concluded 
that Hawking's derivation does not seem to be altered by vacuum fluctuations of linearized 
quantum gravity. Parentani lH and Barabes et al jlH have discussed the possibility that 
quantum fluctuations could be the source of the non-linearity needed for the mode regener- 
ation picture. 

Most of the papers cited in the previous paragraph deal with active fluctuations, the 
spacetime geometry fluctuations arising from quantization of the dynamical degrees of free- 
dom of gravity itself. Another source of spacetime fluctuations are the quantum fluctuations 
of matter field stress tensors, which cause passive fluctuations. There has been an extensive 
discussion of both types of fluctuations in recent years in various contexts, including both 



black hole spacetimes 
geometry fluctuations [20|, [21|, |22j, [23|, 




19| and more general spacetime 

12J,|25| 

In the present paper, we will examine some examples of both active and passive fluc- 
tuations. However, our main interest will be in the possibility of enhancing the geometry 
fluctuations above the vacuum level by use of gravitons or matter fields in squeezed vacuum 
states. We consider a Schwarzschild black hole formed by gravitational collapse, and then 
suppose that wavepackets of gravitons or matter fields in squeezed states are sent across 
the horizon. This will cause greater geometry fluctuations than would occur in the vacuum 
states of these fields. The key question which we wish to address is whether these fluctua- 
tions significantly alter the outgoing modes which will carry the thermal radiation to distant 
observers. The technique which we employ to study the effects of geometry fluctuations is 
based upon the geodesic deviation equation. This allows a gauge invariant treatment using 
the Riemann tensor correlation function. 

In Sect, mi we review Hawking's derivation of black hole evaporation, and discuss the 
transplanckian issue. In Sect. IHIl we develop some of the formalism of fluctuations of a 
geodesic separation vector which will be used in subsequent sections. We next turn to the 
case of active fluctuations. Before considering gravitons, we first investigate a simplified 
model of "scalar gravitons" in Sect. IIVI This model reproduced the essential physics of the 
effects of gravitons, but with reduced complexity. The case of gravitons, quantized linear 
perturbations of the Schwarzschild geometry, is treated in Sect. |V] We next turn to passive 
fluctuation effects in Sect. PVIl where stress tensor fluctuations of a scalar field are treated. 
We give a unified analysis of the results of all three models in Sect. IVHl and offer our 
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collapsing matter 



FIG. 1: The spacetime of a black hole formed by gravitational collapse is illustrated. The shaded 
region is the interior of the collapsing body, the r = line on the left is the worldline of the center 
of this body, the r = line at the top of the diagram is the curvature singularity, and is the 
future event horizon. An ingoing light ray with v < vq from J?~ passes through the body and 
escapes to as a u = constant light ray. Ingoing rays with v > v q do not escape and eventually 
reach the singularity. 

conclusions in Sect. IVIIII 



II. DERIVATION OF THE HAWKING EFFECT 

In this section, we will briefly review Hawking's derivation [lj of black hole evaporation. 
The basic idea is to consider the spacetime of a black hole formed by gravitational collapse, 
as illustrated in Figure [TJ Here v — t + r*, and u = t — r* are respectively the ingoing and 
outgoing Eddington-Finkelstein coordinates, also referred to as the advanced and retarded 
times, and r* = r + 2M ln(^g — 1) is the tortoise coordinate. A quantum field propagating in 
this spacetime is assumed to be in the in- vacuum state, that is, containing no particles before 
the collapse. In the case of a massless field, a purely positive frequency mode proportional to 
e -iujv i eaves propagates through the collapsing body, and reaches J? + after undergoing 
a large redshift in the region outside of the collapsing matter. At J^ + , the mode is now 
a mixture of positive and negative frequency parts, signalling quantum particle creation. 
Of special interest are the modes which leave J"~ just before the formation of the horizon, 
which is the v = Vq ray. These modes give the dominant contribution to the outgoing flux 
at times long after the black hole has formed. After passing through the collapsing body, 
they are u = constant rays, where 

„ = _4Mln(^), (1) 

where M is the black hole's mass, and C is a constant. The logarithmic dependence leads 
to a Planckian spectrum of created particles. It also leads to the "transplanckian issue", 



the enormous frequency which the dominant modes must have when they leave J>~ . The 
typical frequency of the radiated particles reaching J :+ midway through the evaporation 
process is of order 1/M, but the typical frequency of these modes at is of order 

lo « M _1 e (M/Mpi)2 , (2) 

where Mpi is the Planck mass. Another way to state this is to note that the characteristic 
value of u for these modes is of order 

u c *m(-J . (3) 

A geodesic observer who falls from rest at large distance from the black hole will pass from 
u = u c to the horizon at u = oo in a proper time of 

5t^M e ~ Uc/4M w M e - M2/m i . (4) 

which is far smaller than the Planck time. In this sense, the outgoing modes are much less 
than a Planck length outside the horizon. 

If spacetime geometry fluctuations cause these outgoing modes either to be ejected pre- 
maturely, or to fall into the singularity, then the outgoing radiation, and possibly the thermal 
character of the black hole, could be greatly altered. This is the question which we wish to 
address in this paper. 



III. FORMALISM 



A. Null Kruskal Coordinates 

Most of this work is done using null Kruskal coordinates, for which the line element is 

ds 2 = -^Ml e ~ r / 2M dUdV + r 2 dQ 2 (5) 
r 

where the coordinates (U, V) are defined by 

U= _ e -u/4M and y = e v/AM (g) 

and describe surfaces of constant phase, equivalently the path of light rays, in a Schwarzschild 
space-time. Kruskal coordinates are advantageous because, unlike Schwarzschild coordi- 
nates, they do not suffer a coordinate singularity at the horizon. 

In null Kruskal coordinates, V is an affine parameter on M' + and is approximately an 
affine parameter on outgoing null geodesies very near the horizon, but only on that portion 
of the geodesic for which UV <C 1, or equivalently near the r = 2M surface. However, 
outgoing null geodesies spend a long affine time near r = 2M before finally escaping to 
infinity, so this should be a good approximation for a large range of V. 

Furthermore, on the past horizon of an eternal black hole, Jrf?~, the Kruskal coordinate U 
is an affine parameter for ingoing null geodesies. (See, for example, Eq. 12.5.10 in Ref. [HJ.) 
Working in Kruskal coordinates and using U and V as affine parameters near the past and 
future horizons, 1^ = (0,1,0,0) and = (1,0,0,0) are tangent to, respectively, outgoing 
and ingoing null geodesies near the horizon. 
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collapsing matter 



FIG. 2: Schwarzschild black hole space-time formed via gravitational collapse. The separation 
vector characterizing geodesic deviation of the horizon, and a nearby outgoing null geodesic, 

7, is initially n u = Uq but then evolves as n v = 5 s n^. The vector is tangent to the horizon. 
Perturbing fields fall into the black hole well after the horizon forms. 

B. Geodesic Deviation 

The derivation of the Hawking effect outlined in Sect. [Ill requires propagation of a field 
from backwards along a geodesic, through the collapsing body, and out to . These 
tracked geodesies lie very close to the horizon, and are separated from the horizon by some 
separation vector n M as in Fig. El In this figure, the geodesic (labeled 7) appears to be a 
straight line at a fixed distance from the horizon so that is constant. This is not quite 
true; a particle following 7 is eventually separated from the horizon by an infinite physical 
distance. Tracking the evolution of the separation vector from some initial point out to 
requires integration of the geodesic deviation equation from some initial point out to J^ + . 
One finds that the U component of n M is constant while the V component is not. Hawking 
is actually considering only the fixed U component as it is only n u that is relevant for his 
derivation of Eq. ([I]). 

Consider parametrized by U — and affine parameter Ai = V, and a nearby 

outgoing null geodesic just outside the horizon parametrized by U = and affine parameter 
A 2 ~ V for some small no <C 1. Let connect points of equal affine parameter on J^f + and 
the geodesic with U = n . Parameterizing the geodesies such that initially Ai = A 2 = V , the 
separation vector is initially Uq = (n , 0,0,0) at some point V = V . In a flat space-time the 
separation vector would not change and n fl = Uq everywhere along the geodesic. By applying 
the geodesic deviation equation in Kruskal coordinates one can find the subsequent evolution 
of the separation vector = rtQ + 5 s n^, where 5 s n^ represents the kinematic evolution of n M 
due to the classical background. If in addition, there is a perturbation of the background, 
then n M = n,Q + 5 s n^ + <5 p n M , where 5 p n^ represents the dynamical response of to the 
perturbation. The bar on n M is used to differentiate between the background-only space- 
time separation vector and the background plus perturbation space-time separation vector 
for this discussion. 



1. Background Space-Time 



Consider first the unperturbed Schwarzschild space-time of Fig. [2j Let = (0, 1, 0, 0) be 
tangent to the horizon, and let = Uq + 5 s n^ denote the separation between the horizon 
and a nearby outgoing null geodesic with initial separation Uq at V = Vq such that rift is null 
with ri^l^ = 1. It suffices to choose nft = (n , 0,0,0) with n = (guv)~ 1 \r=2M = 1/(8M 2 ). 
The evolution of the vector characterizes the geodesic deviation of these outgoing rays, 
and obeys the set of differential equations 

^ = ir^m. (7) 

Since £ M = (d/dVY, the geodesic deviation is 

D 2 n a ^ D 2 (5 s n a ) 
dV 2 ~ dV 2 
The only non-zero component of R a VVu is 



K-wfFn" . (8) 



UVr 3 

which near the horizon reduces to 

R v vvu « -2c- 1 . (10) 

Also near the horizon, the covariant derivative with respect to V on the left hand side of 
Eq. (jSJ) reduces to an ordinary derivative, which may be seen by direct calculation. The 
second covariant derivative is 

^ = ffn" iff7 + T^P£ 5 n x + 2Y^H s n\ 1 + T^ sx P£ s n\ (11) 

Each term with a Christoffel symbol is of the form Vy . A straightforward computation of 
the Christoffel symbols in null Kruskal coordinates reveals that T v — > as U — > for all a 
and fi. The only non-trivial equation is then 

d 2 (5 s n v ) n , 

= -2e~ 1 -o (12) 

which may be integrated from the initial point V = Vq, resulting in 

n" = n (1, -e~\V - V Q f) or (n^) 2 = 4 (m^" 1 ) 2 {V - V ) 4 . (13) 



2. Perturbed Space-Time 

We wish to let the space-time fluctuate in some way and describe what happens to the 
outgoing null geodesies. Consider an ensemble of Schwarzschild space-times such that the 
average apparent horizon is defined by r = 2M for some chosen value of M. Consider 
an outgoing geodesic just outside the average apparent horizon. For each space-time in the 
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ensemble, define the separation vector n M such that in the average space-time = n M where 
n M is the classical solution for the Schwarzschild space-time of mass M. 

Let g^y = + h^ v be the perturbed space-time, where g^ v is the unperturbed 
(Schwarzschild) background space-time and h^ u is the perturbation. The background metric 
is used to raise and lower indices. The Riemann tensor is calculated from 

td& pa _ po _j_ pet -per -pa per / -i a\ 

where 

f % = ^p a<i (9^,u + gpu,„ - g^,p) (15) 

are the connection coefficients in the perturbed space-time. These may be expanded to 
r™ = r" + 5r™ , where T™ v are the connection coefficients of the background, and 5T^ U is 
due to the perturbation. This yields 

to first order in the metric perturbation, where similarly B? 1 ^ and bR a p denote the 
background and perturbation contributions to the Riemann tensor. Let a semicolon denote 
covariant differentiation with respect to the background One may then verify that 

$R a ^v = ^ r ^; M - $ T ^;u , (17) 

where 

Let £ M be fixed as tangent to the apparent horizon in the average (or background) space- 
time. Let = + 5 s n^ + 5 p n^ denote the separation vector, where Uq is the same initial 
separation as before, <5 s n M is defined to satisfy the background equation as above, and 5 p n^ 
encodes the dynamical response of n M to the perturbation. Letting D be the covariant 
derivative with respect to the perturbed space-time, the geodesic deviation is 

D*n« _ D\5 s n«) D\5 p n") _ p 

dv* ~ dv* + dv* ~ u ^ 1 n 

= (ifV + ^V) W « + W + 5 p rf) . (19) 

Consider the left hand side of this equation. The background terms involving F v still 
vanish near the horizon. There are, however, terms involving 5Y Vfl which do not vanish on 
the horizon and which enter to first order in the metric perturbation. In particular, the first 
term on the left hand side is 

£>2 y ) = m<^), CT 7 + ST^P£ s (5 s n x ) + 25Y^i\5 s n x ) n + 5T^5T^P£ s (5 s n x ) , 

(20) 

which contributes two terms to first order in the metric perturbation. Similarly for the 
second term on the left hand side: 



D 2 (5 r& 

(21) 
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However, 5 p n x is already first order in the metric perturbation, so terms involving both S p n x 
and are second order. Therefore to first order in the metric perturbation, the left hand 
side of Eq. ([19]) is 



AM 



= f^T + f^, + + 25r ^ {S ,n% , (22) 

and the right hand side of Eq. ( fT9l) is 

("o + <W) + BTtfuPPiSpn?) + SR\ P J^ K + 5 s n») . (23) 

The first term is just the result obtained for the background. Since 5 s n^ is by definition the 
solution to 

= R\,J^ K + <U") , (24) 

these terms cancel. Furthermore, since 8 s n^ has only a V component and 5i?°yyy = by 
the antisymmetry of the Riemann tensor on the last two indices, then 8R a ^ u ^^{5 s n v ) = 
and Eq. (TT9~j) becomes 

= R\ p J^{5X) + 5R\ pv Pfin» - ST^P£ s (6 s n x ) - 26T^P£ 5 (5 s n x ) >T (25) 

The solution to this equation involves integrating twice over V from some initial point 
Vo. It has already been found in Eq. ([TBI that 5 s n x oc (V — Vo) 2 . We suppose that 5 p n x 
may also be expanded in powers of (V — V ) and then approach the solution to Eq. fl25|) 
iteratively. Suppose 5 p n x = 5 p ^n x + 5 Py 2n x in powers of (V — Vo). Since the constant term is 
already accounted for in n x , then 5 p ^n x must be 0(V — Vo) or smaller. The double integral 
over the 5 p n x and 5 s n x terms results in terms of higher order in (V — Vo), thus in the first 
iteration, 5 Pt in x will be the solution to 

= SR\?J^<, (26) 

which turns out to be proportional to (V — Vo) 2 . The next iteration is the solution to 

d 2 (6 Pt2 n a ) 



dy2 ~ R\ Pv ^{5 Pjl n») - 5Y^£\5 s n x ) - 25Y^£\5 s n x ) , r (27) 

Each term on the right hand side is proportional to (V — Vo) 2 , and the solution is simply a 
double integral over V, giving a result proportional to (V — Vo) 4 . Consequently, to lowest 
order in powers of (V — Vo), the solution to n M for the perturbed space-time is 



pV j-W 

n° + 5 s n a + dW dV5R X tfjn p r%. (28) 



To proceed further requires a model for fluctuations to be specified, and here three dif- 
ferent models will be considered in turn: 

• A scalar graviton model where the metric is perturbed by a term proportional to the 
product of a scalar field with the background metric, i.e. g^ u = (1 + $)^. 
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• An ingoing gravitational wave actively perturbs the horizon in a linearized theory of 
gravity. The perturbation to the Riemann tensor arises from the perturbation to the 
metric. 

• An ingoing scalar field provides a passive perturbation to the horizon. The perturba- 
tion to the Riemann tensor is the Ricci tensor contribution that arises from the stress 
tensor of the scalar field. 

In all three cases, the ingoing field is taken to occupy a squeezed quantum state, |a, £)■ 
In particular, the expectation value will be evaluated with respect to a multimode squeezed 
vacuum state, |0, () = YliLzo ^(Ci)|0), which is further described in AppendixO The excited 
modes will be taken to be wavepackets which are sent into the black hole after the collapse, 
as illustrated by the ingoing arrow in Fig. [2j 

It should be noted that the introduction of quantum fluctuations into a black hole space- 
time entails a significant conceptual shift. Classical perturbations will shift the location of 
the horizon, but do not change the fact that there is a precisely defined horizon. Of course, 
the true event horizon in a classical spacetime, the light ray which barely fails to escape to 
J^ + , can only be known when the complete history of the spacetime is known. Quantum 
fluctuations introduce an additional ambiguity, whereby the precise event horizon can never 
be known. 



3. Fluctuations 

Quantizing the ingoing perturbation field, 5 p n^ becomes a quantum operator. To con- 
struct the operator S p n^, consider Sph^ as the solution to 

^ = Sfr^Prt. (29) 

Suppose we characterize fluctuations in the separation vector by the quantity 

(n^fifi) - (n^)(n u ) = (5 p n >t 8 p n lt ) - (c^XVv). (30) 

Due to the peculiarities of null Kruskal coordinates, this quantity is not a good comparator 
for all three models. In particular, it is identically zero for the "scalar graviton" model, of 
order (V — Vq) 4 for the graviton model, and of order (V — Vq) 2 for the scalar field model. In 
order to compare the results of the three models with each other and with the background 
result, which is of order (V — V ) 2 , it is advantageous to consider instead the variance 

A(r^rv) 2 = ((n^x)n^x)) (n^x')n^x'))) - ((n^x)n^x)))((n^x')n fl (x'))). (31) 

It is straightforward to show that 

A(n%) 2 = 4 «(xK(x') + n%(x)5,n"(x f ) + 5 s n»{x)5 s n u (x')) 

x [(Sph^Sphvix')) - (5 p n^x))(5 p n„(x'))] + 2 (nft(x) + 5 s ri\x)) 
x [{5 v hp{x)8 p h v (x')S p n u (x')) - (Sph^x)) (5 p h v (x')5 p h u (x'))} 
+ (5 p h^(x)5 p h fl (x)5 p h' / (x')5 p h u (x')) - (5 p h^(x)5 p h^(x))(5 p h u (x')5 p n u (x')) . (32) 
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As will be demonstrated, each 8 p h v {x) oc (V — Vo) 2 ; additionally, S s n v (x) oc (V — Vo) 2 . 
Therefore, to lowest order in (V — Vo) we have 



A(n M n M ) 2 = A[{(n^(x)8 p n Pj (x))(n v (x')8 p h v (x'))) -{n^{x)5 p h ll {x)) (n u (x')8 p h u (x'))} . (33) 



This then is the primary quantity of interest to calculate for the three fluctuation models. 



IV. SCALAR GRAVITON MODEL 

Turn now to the scalar graviton model, where the perturbation is simply a scalar field, 
$, (a dilaton) multiplying the background metric, = (1 + $) g^. Here g^ is the unper- 
turbed (Schwarzschild) background space-time metric and we may define g^ u = g^ u (l — $) 
such that to first order in $, g^ a g^p = 8 a f3 . The scalar field $ will be a free quantum scalar 
field (multiplied by ip) and will be defined in the following section. One may object that 
this model is simply a conformal transformation of the space-time under which the light cone 
structure remains invariant. However, while it is true that the light cone is invariant under 
a conformal transformation, the geodesic deviation is affected. This is because the Riemann 
tensor involves derivatives of the conformal factor, so that the Riemann tensor of the trans- 
formed space-time is not equal to a simple conformal transformation of the Riemann tensor. 
For example, Robertson- Walker space-time is conformally flat but has non-trivial geodesic 
deviation. This model is useful as a simplified model which reproduces the essential features 
of the more complicated graviton model of Sect. [V] 



A. Normalized Wave Packets 

Let a scalar field propagate from r* = oo, through the potential barrier of the black 
hole, to the horizon at r* = — oo. The wave function must satisfy the wave equation in the 
Schwarzschild geometry. Following Hawking Fourier decompose solutions of the wave 
equation with respect to advanced or retarded time, use continuum normalization, and 
expand in spherical harmonics. Using ingoing Eddington-Finkelstein coordinates, a single 
ingoing mode is 



The integer j controls where in frequency space the wave packet is peaked, while Ej controls 
the width of the wavepacket and has units of frequency. The integer n describes which 
wave packet is under consideration. This construction allows for wave packets to be sent 
in at regular intervals of 2ir/ej with various frequencies. Thus ipj n is the n th wave packet 
sent in with component frequencies ranging from jej to (j + 1)ej. The function F w {r) is 
in general a complex function which depends in some complicated way on the geometry of 
the space-time. For sharply peaked wave packets, F^ir) is of order unity at infinity, and 
essentially reduces to a transmission coefficient near the horizon. In Appendix |A], it is shown 




(34) 



from which ingoing wave packets may be constructed as 




(35) 
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that these wavepackets are properly normalized. The quantized scalar field is constructed 
from the wavepackets defined above as 

$ = Yl (f^a jn + V&ojn) • ( 36 ) 

jn 

In general it suffices to consider a single ingoing wavepacket, thus in what follows the index 
n will usually be suppressed and assumed fixed. 



B. Fluctuations 

Since the metric of the full space-time obeys the same symmetries as the background 
space-time, it is straightforward to calculate the Riemann tensor exactly from Eq. ( fT4l) . In 
particular, to first order in $ and its derivatives, one finds the relevant quantity [see Eq. ( fTBT) ] 

5R V VVU = ® yu (37) 
Applying Eq. (1281) one finds, to leading order in (V — Vq), 

n^x)5 p n^x) = g uv n 2 dW dV 5R V VVU = g uv n 2 dW dV$ >vu . (38) 

JVo Jvo Jv JVo 

Recall that guy = 8M 2 near r = 2M. Clearly, ($) = even when evaluating the expectation 
value with respect to a squeezed vacuum state; therefore the variance, Eq. ( |33|) . becomes 

A(n%) 2 = 4(«(x)Vv(x)) {n»{x')5 p n v (x'))). (39) 

Expanding in terms of the mode functions, one finds 

-V pW pV pW 

dW dV dW I 

'V JVo JVo JVo 



pV pVV pV pW 

A(n^) 2 = 4(g uv nl) 2 Y] / dW dV dW / dV 

j:k Jv JVo JVo JVo 

x \$j,vu{x) ipk,vu'(x') (C, 0\aja k \0, C) + ipj,vu(x) V k yu'( x ') (C, 0\djd{\0, () 

+r h vu^)^yiu>(x') (C,0|ata fc |0,C) +V'*v C /(^) ^tyu'W) (C,0|atat|0, C>) (40) 



where we choose to evaluate the expectation value with respect to a multimode squeezed 
vacuum state |0, £) = YliL Z0 S(Ci)\®)- Using the results of Appendix ICl this becomes 



pV pW pV pW 

A(n%f = 4( 9w ^) 2 V / dW dV dW / dV 

, , JVn JVn JVn JVn 



j,k JV ° JV ° JV ° JV ° 

x 5 jk {i) j y U (x)i) k yi U i(x') (1 + Q ZoZ1 (j) (cosh Pj - 1) (-0 2o2l (k) sinh p k ) 
+ ipjyvix) Tp* k yu>(x') (1 + 2o2l (j)(coshp.,- - 1) (1 + 6 2()2l (£;)(coshp fc - 1) 
+ ^j,vu( x ) ^k,vu'(x') (-0 2o2l (j) sinhpj) (-6 2o2l (/c) sinhp fc ) 
+^j,vu( x ) ^*k,vu'( x ') (-@202i U) smh pj) (1 + 9 202l (k)(coshp k - 1)} , (41) 



11 



with the integer step function 



®Z Z! (j) 



1, z <j< Zi, 
0, otherwise. 



(42) 



Notice that the factor multiplying ^jyu{x)'4 ) kV'U'( x ') contains a 5jk which makes 
A(n /t n M ) 2 divergent. We therefore take renormalization to correspond to restricting the 
sum over modes to those occupying an excited squeezed state mode. For the current 
model and the graviton model, this restriction corresponds to normal ordering. In gen- 
eral, there will also be a vacuum contribution which is being ignored here. This should be a 
good approximation for highly excited states, that is, the limit of large squeeze parameter, 
coshpj pa sinh pj pa e Pj /2. Taking this limit and restricting the sum to zq < j,k < Z\, the 
integer step function Q Zf)Zl = 1 and this simplifies to 



zi „V pW 

A(n^) 2 = {g uv nl) 2 V e« dW dV (^ yu (x) - ^* yu (x)) 

j=z J v° J Vo 



The V integral is trivial and we have 



(43) 



3=zo 



U hU\ v=Vo ^3,U\ V =V 



(44) 



Consider ipjjj(x). Near the horizon, F u (r) reduces to a transmission coefficient which de- 
pends only on uo and £. One may then find 



(3'+l)ej 



duj 



* im.fi d J _ l/ \ T 7-1 — AiMu) 



tF(uj)V 



UJ 



(45) 



where 8j = 2iin/ej. We recall that (1 — — ) ~ —UVe 1 , and use the definition of V to write 



Ve -iu>v = v i-4iMu>_ Next) one findg 



V 



dW^x 



\V=W 



duo 



Y lm F(u)e 



! u ./, . " 2.\/ x /2 
x [(2 - AiMuj)' 1 (V 



71UJE 



2-4iMu) 



J 

v 



(46) 



Expanding the bracketed terms in powers of (V — Vq), this becomes 



/; 

Jv 



dWipj^ix 



v=w 
v=v 



'J'i 



(47) 



Here Vq is the Eddington-Finkelstein coordinate corresponding to Vo, i.e. Vq = e v °^ M . Use 
Eq. f|4"T|) in Eq. (144")) and let m = 0, then Y^ = and the spherical harmonics factor out. 
Using this and the results for the classical deviation, Eq. (fT3|) . and inserting the appropriate 
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powers of the Planck length, £p, the fractional fluctuations are found to be 



(n^n u ) 2 ^ 
p i=zo 



du 



eo<-p 



8M v /2vfcJi~ 



x [(F(u)e- Mvo+5 ^ - F*(uo)e iu{vo+ ^) - AiMu (F(cj)e- luj( - V0+5 ^ + F*(u)e iuj{vo+5 ^)] 



(48) 

Technically, expanding Eq. (T4"6T) in powers of (V — Vo) is really an expansion in powers of 
2Muj{V — Vo). In considering the result for the fractional fluctuations, therefore, one must 
bear in mind that we are considering the limit where (V — Vq) < 1 and also to < (2M)~ l . 
In the limit u — > 0, F(u) — > B^{2iuMY +1 where \B(\ is of order 1. In this case one finds 



A(n'X) 2 _ \B e \ 2 Y 2 A e 2 ^£ 2 P (2M) 2e 
(n»n^) 2 8 ^ 



TTEj 



du)uj e+1 / 2 cosu;(fo + Sj) 



(49) 



Near uj ~ (2M) 1 , the transmission coefficient is of order 1. It follows that if the wavepacket 
is sharply peaked near uj « (2M) _1 



A(n^n M ) 2 _ Y 2 y, e 2 ^t\ 



16 ^ MTie j 

j = ZQ J 



0'+l)ej 



do; sinu;(t>o + Sj) + 2 coscj(t>o + 5j 



(50) 



Let us leave further analysis until Section IVIIII and first consider the other models. 



V. GRAVITON MODEL 

The next model for fluctuations is that of an ingoing gravitational wave occupying a 
squeezed state. In this case the perturbation field mode functions are constructed from the 
allowed classical black hole perturbations. 



A. Even Parity Classical Perturbations 



With Regge and Wheeler paving the way, the subject of classical black hole perturbations 
was thoroughly studied by Vishveshwara, Eddlestein, Zerilli, Price, Tuekolsky, and others 



29l I30L 3ll . I32I [H, [H] . We begin with the original formulation of metric perturbations 



by Regge and Wheeler. These come in two varieties, even and odd parity. In this work, 
we will give an explicit treatment for the even parity case. However, gravitons in odd 
parity wavepackets can be shown to lead to similar conclusions as we will find here. Purely 
even parity waves are physically realizable, being generated, for example, by matter falling 
radially into a black hole [HI]. I n Schwarzschild coordinates and using the Regge- Wheeler 
gauge, the even parity metric perturbation is 



—iu)t 



fiu 



PA cos t 



Hi{r) (1-^)- 



1 H 2 {r) 



r 2 K{r) 



\ 



! sin 2 9K(r) ) 



(51) 
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While this equation pertains to a particular choice of gauge, our results are based on the 
Riemann tensor correlation function and as such are gauge invariant. Zerilli 30|, |31( showed 
that the even parity radial functions H (r), Hi(r), H 2 (r), and K(r) may be related to a new 
radial function, Z(r), that obeys a single Schrodinger-type equation 

^ + - V eff )Z(r) = 0, (52) 

with an effective potential 

(I- 2M\ 2A 2 (A + l)r 3 + 6X 2 Mr 2 + 18AM 2 r + 18M 3 
Veff ~ \ ^~) r3(Ar + 3M)2 ' (53) 

where A = (£ - 1)(£ + 2)/2. 

Equation (l5Tj) may be transformed to null Kruskal coordinates, which we indicate with 

a superscript, \I>^/. Next, construct the wavepackets 

*£U = /. dwiW^e-^tfg?, (54) 



where Sj = 2im/ej and A£ m (aj) is a normalization factor. The normalization of this pertur- 
bation mode requires some amount of work, which is relegated to Appendix [B] The ingoing 
linearly quantized graviton perturbation field is then 

« ) = E (*83^» + HC ) • ( 55 ) 

jn 

The interpretation of the integers j and n is exactly the same as discussed in Sect. IIV Al 
Again, it suffices to consider n fixed, so the index n will in general be suppressed. Using the 
results of Appendix [Bj the properly normalized even parity wavepacket is written 

*tfW=/ ^=e"^^, (56) 



where 

1 = 2ZTI [2 + f(i + m + 3)] - f^lf ■ (57) 

With this normalization, it is understood that the Zerilli radial function takes the asymptotic 
value Z(r* — > oo) = e~ lujr * 



B. Fluctuations 



From Eqs. (fTBl and (I55I) . one finds 

= \^ E /° +1)£j rfw [^m(w)e-^ (tt^ + - <Ew) o i + H.C.] , (58) 

2 3 J ™ 
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and by extension 



M2V = ]T /° +1)£j du [A im (u)e-^SR' Pltv a j + H.C.] . (59) 

Here the hat has been reinserted to clarify that the operators ST^ U and 8R a p may be 
constructed from the classical quantity corresponding to a classical single mode perturbation. 
Consequently, 

S p n» = J2 du} [Atm(w)e '^0.3 + H.C] . (60) 



The second term of Eq. (1331) is zero for squeezed vacuum states. Recognizing that Ag m (uS) 
is real, the variance for the graviton model is 

A(n M n M ) 2 = 4((n£(x)5 p h p (x)) {n v Q {x')8 p n u {x'))) = 

rti+ifo r(k+i)e k 
^Z2z2 duJ du'A im (u)A im (uj') 

| e - iM ,W4) ( n »{x)5 p n,{x)) {n v Q {x')8 p n v {x')) (dfa k ) 
+ e^+^W ( n $(x)8 p n;(x)) K(x%<(*')) (a)a\) 
+ e-^-^ «(x) W(x)) (n»{x')5 p nl{x')) (afa\) 

+ e^~^ (n%(x)6 p n*Jx)) (n u (x')S p n v (x')) (a}a k )\ . (61) 



We choose to evaluate the expectation value with respect to a multimode squeezed vacuum 
state \0,C) = YliLz 'S'(O) 1 0) . As discussed for the scalar graviton, renormalization amounts 
to restricting the sum to those states which lie in the range of squeezing. Together with the 
results of Appendix O and in the limit of large squeeze parameter, p, this gives 
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a(^) 2 = Y, e2pj 

j=ZQ 



duA em (uj) (e-^n%(x)5 p n^x) - e^n^(x)5 p n;(x)) 

3 e 3 



(62) 



It remains to calculate n^pfi^ due to a single mode classical perturbation. Unfortunately 
this is difficult to do analytically, but is possible with the use of a computer algebra program. 
To further simplify the problem we restrict attention to perturbations of purely even parity 
with i = 2. To calculate 5 p n fl near the horizon, one must first calculate SR 01 ^^ near the 
horizon in null Kruskal coordinates. There are several possible routes toward obtaining this 
information. The most straightforward might appear to be to first transform the metric 
perturbation to Kruskal coordinates and then proceed from Eq. (IT7|) . The difficulty with 
this method arises when one tries to take the limit of SR a g near the horizon. In null 
Kruskal coordinates, the approach to the horizon is along a line of constant V rather than 
a line of constant t. The equations, however, still involve r, implicitly defined in terms of U 
and V, so taking the limit is not a well defined operation. 

It is instead simpler to begin with the metric perturbation in Schwarzschild coordinates 
and compute SR ^^ via Eqs. (fTTI) and (|T8l) for a classical single mode ingoing perturbation. 
Expressing the radial functions H(r), Hi(r), and K(r) in terms of the Zerilli function, Z(r), 



15 



one then uses Z(r) to first order in (r — 2M) to expand 5R a 'a in powers of (r — 2M) near the 
horizon. It turns out to be necessary to expand the Zerilli function to at least 0(r — 2M) 
to ensure that the metric perturbation remains finite on the horizon. Since the initial 
calculation is done in Schwarzschild coordinates, the limit is well defined as r —>■ 2M along 
a line of constant t. Once SR ^^ is known near the horizon, it may then be transformed to 
null Kruskal coordinates. 

Using this procedure to solve for 5 p n^ as expressed in Eq. fl28|) . one may go on to find 



n^(x)8 p n lx (x) = -{guve n )(V - V ) 



T(uj)e- iuJV ° 



(11800i 



324M(z + AMu) 
+34321Mcj - 39783iMV + 14348MV) (1 + 3 cos 26) . (63) 

Inserting the appropriate powers of the Planck length, the fractional fluctuations are then 



i=zo 



(i+ife £ P (1 + 3 cos 26) \T(u)e~ iul{vo+5 ^ 



324M\/ire j Luj 



(i + AMu) 

(11800* + 34321Mcj - 39783zMV + 14348MV) - H.C. 



• (64) 



Once again, one must consider that the expansion in (V — Vq) is really an expansion in 
2Muj{V — Vq). In the limit u — > 0, the transmission coefficient is T(u) ~ Cg(2iuM) +1 , 
where similarly to the scalar case, \Ce\ 2 ~ 1. Specializing to I = 2, T(u) ~ — iC2{2M) 3 uo 3 



and L — ^P. Then 

5 



A(n^) 2 _ 5 2 / 1475(1 + 3 cos 20) \ 2 ^ e 2 «^(2M) 4 
\ C *\ { si ) 2. ^ 

' ]=ZQ 



(n^n^) 2 ~ 2 



7T£ 1 



(J'+I)ej 



1 2 



duuj 5 / 2 cosco>(t> + 5j) 



(65) 



When u pa (2M) 1 , the transmission coefficient is of order 1. It follows that if the 
wavepacket is sharply peaked near u pa (2M)" 1 , then 



A(n^n^) 2 (l + 3cos2#) ; 
(n^) 2 
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^ e 2 ^£ 2 P 1 



7T£j M 



0'+l)ej 



rfo; 21.1 sin o;(fo + 5j) + 54.9coso;(vo + 5,-) 



(66) 



The physical content of this expression will be explored in Sect. IVIIIl 



VI. PASSIVE FLUCTUATION MODEL 



This model differs from the scalar graviton and graviton scenarios in that rather than a 
quantization of the dynamical degrees of freedom of the gravitational field, the space-time 
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geometry fluctuations arise passively through fluctuations in the stress tensor of a quantized 
scalar field. The ingoing scalar field is constructed in the same manner as presented with 
the scalar graviton model. The quantity of interest is the variance of the squared length of 
the separation vector, A(n M n^) 2 , given by Eq. (1331) . Unlike the two previous models, the 
second term of this equation is not zero for the present model. 

The operator nature of &R? av p is due to the stress tensor of the scalar field. One can use 
(see e.g. Ref. |2fj|, Eq. 3.2.28) 



R 



afifiu 



(n-l)(n-2) 



with 



to find 



R 



^ ( - -g^T ) , and R = -8<kT 



5R, 



Ol/3fMU 



8tt 



(67) 



(68) 



(69) 



Here the perturbation of the Weyl tensor vanishes, 5C a p^ u = 0, and the number of space-time 
dimensions is n = 4. The stress tensor for a scalar field and its trace are 



T, 



$>^v) - 2^ CTP$ ;(^;p)> and T = -9 ap ®;(<r®; P )- 



It follows that 

oR a j3/j,u 



•>7T 



where it is to be understood that the antisymmetrization proceeds first, i.e. 



To leading order in (V — Vq), 



™o<Vv = 9uv n 2 



dW 



Vo 



dV5R l 



vvu- 



(70) 



(71) 



(72) 



(73) 



The variance is therefore proportional to an integral of a component of the Riemann tensor 
correlation function 



pV pW pV pW 

A(n%) 2 = A{g uv nlf / dW I dV dW / dV (C 

Jv Jv Jv Jv 



i iriV V 



VVU VVU 



(x,x')) (74) 



where 



(C v wu V wuM) = (5R v vvu (x)5R v vvu (x')) - (5R v vvu (x))(SR v vvu (x')). 
From Eq. (1711) . the Riemann tensor component of interest is 



(75) 



5R 



v 



8tt 



VVU 



($ ; y$ ;C 7 + ®;U®;v) ~ ~^UV [<T (®;0$;0 + + g W + A 



(76) 
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For simplicity we restrict to the case £ = 0, so the angular derivatives are zero. Expanding 
$ in terms of its mode functions gives 



V vvu — i 



+ {^lu^k,v + ^j,v^k,u) a}a k + (^,^*,v + ^j,v^l,u) aJ-4J ■ ( 77 ) 
Simplify the notation by writing 

(J-R^yy^ = ^ j^(x)aja fc + E^a^ + B*(x)a} j a k + A*(x)a^a\ (78) 



with v4*(x), B(x), and _B*(a;) defined by comparison with Eq. (1771) . 

To find the expectation value with respect to the multimode squeezed state |0,£) = 
YliLz ^(Cj)|0), use the results of Appendix O We again take renormalization to correspond 
to restricting the sum over modes to those occupying a squeezed state. In the limit of large 
squeeze parameter, p, one finds 

(£ 0| : SR v vvu (x) :\QX) = -—J2 ~^ 5 ^ ^ + A ^ x ) ~ B ^ ~ ■ ( 79 ) 

j,k=z 

and 



j,k=zo r,s=zo 

x [(A(x) + A*(ar) - - B*(x)) (A(x') + A*{x') - B{x') - B*{x'))\ . (80) 

In general, quartic operator products can be expanded into a sum of a fully normal 
ordered part, a cross term, and a vacuum part. (See, for example, Ref. [12].) Our procedure 
of restricting the sum to those modes which lie in the range of squeezing is, in the limit of 
large squeeze parameter, equivalent to retaining only the fully normal ordered part. To see 
this, consider the difference between one of the terms above and its normal ordered version. 
As a concrete example, consider 

(C, 0|a i aj b a r a s |0, - ((C, 0\alaja r a s \0, () + 5 jk ((, 0|a r a s |0, ()). (81) 



From equations (1C15il) . (1C15ml) . and (1C12aj) it is clear that this is proportional to the 



subdominant term e 2p (compared to e 4p ). Thus, in the limit of large squeeze parameter, 
restricting the sum over modes to those which lie in the range of squeezing corresponds to 
taking the fully normal ordered part. Furthermore, the cross and vacuum terms which have 
been neglected are sub-dominant in this limit. 

The 5jkS rs term is the same as (: 5R v vvu (x) :)(: 5R v vvu (x f ) :), which cancels to leave 



(C V vvu V vvu (x,x')} = ^ (y) E E eftenW&rSk. + 6 JS 6 kr ) 

^ ' j,k=zo r,s=zo 

x [(A(x) + A*(x) - B(x) - B*(x)) (A(x') + A*(x') - B(x') - B*(x'))\ . (82) 
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Using the definitions of A(x) and B(x), noting that the sums extend over the same range, 
and using the fact that the Kronecker deltas act symmetrically, one may show that 



2 21 



(c v wu v wuM) = 4 hH E e2pJe2pk \( l ^(*)),v QnM*)\u 



j,k=z 



X 



(Im^(x')) v , (Im^jfc(x')) v , ■ (83) 



Since (C v vvu v vvu (x, x')) is a product of a function of x with a function of x', the integral 
over x and rr' of the product is the product of the integrals such that 



dW dV dW / dV (C 

JVo Jv Jv Jv 



V V 

vvu vvu 



'x,x')) 
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e 2 pj e 2 pk 



j,k=z 



Vo 



W 



dW / dV(lm^(x)) y (lmM^)),u 



Vo 



(84) 



Near the horizon, the derivatives of ip with respect to U and V are straightforwardly calcu- 
lated. One subsequently finds 



(J'+I)ej 



da; 



(fc+l)e fc 
ke k 



X 



F(a;)F(a;')e" iMjW<5fc) e~™ (w+w,) 

- F(a;)F V)e~ iMj ^' 5fc) e~" ,(t ^' ) + if.C.] . (85) 



The integration over V of (lm(ipj)) v (Im(^)) v reduces to calculating 



fV t-W pV pW i 

dW dVe^ u ^= dW dV V T4tM ^ Tuj,) ^ -e^ v "^'\V - V ) 2 , (86) 

Jv JVo Jv JVo 2 

where the solution has been expanded in powers of (V — Vo). Inserting the appropriate 
powers of the Planck length, it follows that the fractional fluctuations are characterized by 



-v 



-w 



A(n^n M ) 2 



j,k=z 



3ttM 



j du J— (F(u)e- iuj{v " +5 ^ - F*(uj)e ioj( - vo+s ^) 



3 £ j 
■(fc+l)ejt 



X 



/ ad -== (F(Lu')e~ tuj(vo+5k) - F*(uj')e iu>i - m+&k) ) 

Jks k V u ' £ k 



(87) 



Once again, one must consider that the expansion in (V — V ) is really an expansion in 
2Mu(V — Vo). In the limit u — ► 0, the transmission coefficient is again F(u) ~ Be(2iuM) i+1 . 
Specializing to £ = 0, F(lu) ~ 2iB Mu, and the fractional fluctuations are then 



A(n^n M ) 2 
(n^) 2 



t) i b »i 4 E 



: ' e 2pj £ 2 P M 



j,k=z 



US, 



do; o; 3 ^ 2 cosa;(t>o + 5,-) 



x 



7T£ fe 



r(k+l)e k 

j dJ (a;') 1 / 2 cos a/(i> + 5k) 

Jke k 
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Near uj m (2M) 1 , the transmission coefficient is of order 1. It follows that if the wavepacket 
is sharply peaked near uj ~ (2M) _1 , then 



A(n^) 2 _ /4\ 2 ^ e 2 nP P 
(n^n,) 2 ~ V3/ 7T£,-M 

^ 7 j,k=zo J 



do; sinc<j(t> + <5j) 



x 



7T£ fc M 



(fc+l)e fc 



fce fc 



(89) 



VII. DISCUSSION 



A. Summary of Results 

Consider the results for the three different models - scalar graviton, graviton, and stress 
tensor induced fluctuations. There are two limits of interest, uj — ► and uj « (2M) _1 . 



Case 1: uj < 1/M 

In the low frequency limit, — > 0, we found Eqs. (j4~9l . ( 1651) . and ( 1551) for the scalar 
graviton, graviton, and passive fluctuation models, respectively. The solutions to these 
integrals may be expressed in terms of incomplete gamma functions, but it is not necessary to 
invoke the machinery of incomplete gamma functions to get an idea of the general behavior 
of the fluctuations. To simplify the discussion, let us set Sj = 0. This may be assumed 
without loss of generality and is equivalent to assuming n = in Eq. (1351) . In this limit 
we may use the small angle approximation to set cos(ujvo) ~ 1. Furthermore, for a sharply 
peaked wavepacket, we may assume the integrand is approximately constant so that 

J duj uj x cos(ujv ) « uj x Auj. (90) 

Ignoring the numerical factors and recognizing that Sj = Auj, the general behavior of the 
fractional fluctuations is 

Ej ^e 2 ^£ 2 P (2Muj) 2i ujAuj, Scalar Graviton, 
Ej ^e 2 ^£ 2 P {2Muj) A ujAuj, Graviton, (91) 
(Ej \e 2p] ip(2Muj)ujAuj^ , Passive. 



Ain^fiu) 2 



(n^n^) 2 



Case 2: uj (2M)" 1 

For uj ~ (2M) _1 , on the other hand, the fractional fluctuations were given by Eqs. 
(150]) . (!66|) . and (!89|) . Since the wavepacket is assumed to be sharply peaked in uj, then 
Auj = Ej <C uj while j is large. Thus we may use the small angle approximation, cos(EjV ) ~ 1 
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and sin^EjVo) pa EjVQ. Furthermore, as an order of magnitude estimate, we may assume 
sm(j£jVo) is of order 1 so that 

dujsm(uuv ) pa Ej = Auo, (92) 



and similarly 



/ duo cos(cjt> ) pa Ej = Auo. (93) 



The general behavior of the fractional fluctuation in the limit uj pa (2M) 1 is now 



A^n^) 2 J 'Sr^ i Scalar Graviton and Graviton, 

<^1fe^) ! , P-ive. 



Consider for a moment the limit of low squeezing, p — *■ 1. Although our results have been 
derived for p ^> 1, the order of magnitude behavior should still correspond to what we 
have derived since sinh p and cosh p are of order one for small p. In this case the active 
fluctuations behave as £ 2 p Auj/M, which is in agreement with the results of Ford and Svaiter 



11] where they find fluctuations in the proper time of an infalling observer become of order 



one for M ~ l v . 

The scalar graviton and graviton models describe the active fluctuations, while the scalar 
field stress tensor fluctuations are passive fluctuations. Specializing £ = 2 in the scalar 
graviton model, the general behavior is identical to that of the graviton model. In light of 
Eqs. (HHP, ([MD, and (JHZD one sees that Eqs. (JHU), and flMU become 

A(n^n ^ I > — T77T7-4 , Active 




M(Moj) 

(n^n^) 2 | e 2 pe 2 p \F(uj)\ 2 Au 

M 



(95) 
Passive. 



Although we have specialized to £ = for the passive fluctuations, it is not difficult to 
generalize to the more general case of arbitrary £ > 0. For £ > there are some additional 
derivatives of $ with respect to the angular variables, but this will not alter the order of 
magnitude behavior in Eq. ( |95|) . 

There are some generic features common to both active and passive fluctuations. As 
uj — > 0, the fractional fluctuations are suppressed by some power of uj. They are further 
suppressed by a factor of Au, which has been assumed small, and some powers of the Planck 
length. Notice that the passive fluctuations are more heavily suppressed since the passive 
fluctuations are proportional to £ A P whereas the active fluctuations are proportional to £p. 
The suppression in the low frequency limit is expected, since the effective potential barrier 
efficiently prevents low frequency modes from reaching the horizon. However, it seems these 
suppressions may be overcome by arbitrarily increasing the squeeze parameter. 

For uj pa (2M) _1 the similarities between the active and passive fluctuations are even 
more striking. In this case, the fluctuations are again suppressed by the same powers of the 
Planck length and the width of the wavepacket, but are additionally suppressed by the black 
hole mass. One would expect this term to become important once the black hole evaporates 
to the Planck mass. Again, by arbitrarily increasing the squeeze parameter it is possible to 
overcome the various suppressions. 
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B. Semiclassical Restriction on Squeezing 



Since it appears that the fluctuations may become arbitrarily large by unboundedly in- 
creasing p, we should investigate whether there is an upper bound on p. After all, a squeezed 
vacuum state is not necessarily devoid of particle content. From the discussion in Appendix 
IClit follows that for a single mode squeezed vacuum state, |0, £) = £>(£)! 0), the expectation 
value for the number of particles is 

(C,0|A>|0,C) = (01^(0^5(010) = sinh 2 p -> e 2p , p > 1. (96) 

This means that increasing the squeeze parameter increases the mean energy density. If the 
energy density grows too large, then the semiclassical approximation in which backreaction 
is ignored fails. 

Consider the passive fluctuation model. ^From the semiclassical Einstein equation, these 
calculations should remain valid as long as 

« R% u . (97) 

Working in null Kruskal coordinates, a typical component of the background Riemann tensor 
is 

E ^-^r7 J' (98) 

which, since (l — — ) ~ —UV near the horizon, is approximately — 2e _1 . However, 5R V VVU 
is precisely the quantity that was calculated for the passive fluctuation model. Using Eq. 
(1791 it follows that 

(99) 

3 

Consider next the scalar graviton model. In this case SR^^ is linear in the field, and 

therefore (8R a f3nv) = 0. One could look for second order perturbations SR^ oc {h^ u ) 2 that 
would be quadratic in the field. Instead, we consider 

(100) 

to be a good indicator of whether the semiclassical treatment is valid for the perturbation. 
The quantity of interest is 

VVUI ~ V^,VU^,VU 

and the calculations proceed as in Sec. [TV] The result for yu® ,vu) is precisely the same 
as that found for A(n M n At ) 2 /(n /i n At ) 2 in Eq. (HHj) and we therefore have 

(io2) 



For the graviton model, one must again turn to a computer algebra system to calculate 
(5R V VVU SR V VVU ) . The result agrees with Eq. (11021) . In the process of seeking to place an 
upper bound on the results of Eq. (1951) via a restriction on the squeeze parameter, p, a 
restriction on the results of Eq. ( 1951) themselves has been found. It is therefore sufficient to 
proceed with this restriction on the general result. 
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C. Analysis of Results 



For the sake of an order of magnitude estimate of the fractional fluctuations, assume the 
extremal cases where 



E 



e 2 Pil 2 P \F(cu)\ 2 Auj 



M 



1 (103) 



for the active fluctuations, and 

e 2p i£ 2 P \F(uo)\ 2 Auj 



E 



1 (104) 



M(Mu) 

for the passive fluctuations. Using this restriction, the fractional fluctuations become 



A(n^) 2 
{n^n ^ 2 ~ 



< 



1, Active, 



(105) 



^' 1, Passive. 



At first glance, it appeared that the fractional fluctuations described in Eq. f )95|) could be 
made arbitrarily large by increasing the squeeze parameter, a quite surprising result. Since 
n M n M reflects the change in geodesic deviation from the Schwarzschild background, restricting 
to perturbations that are in some sense small would lead one to expect that n^n^ should 
also be small and not deviate very much from the background value, n^n^ (equivalently 
the expectation value of the fluctuating quantity). Indeed, after imposing restrictions on 
the amount of allowed squeezing by requiring the induced curvature to be small compared 
to the background, we find the fluctuations to be no more than of order one. Nonetheless, 
fractional fluctuations of order unity in A(n M n M ) 2 / (n M n M ) 2 have the potential to dramatically 
alter the outgoing radiation. 



D. Implications for Hawking Radiation 

Do the space-time fluctuations implied by Eq. f ll05p have any significant effect on Hawking 
radiation? Fluctuations of the vector separating the horizon from a nearby outgoing null 
geodesic have been studied here precisely because of this vector's importance to the Hawking 
effect derivation. However, it is the U component of the separation vector that is crucial to 
Hawking's derivation. In fact, 5 p n u = by the symmetry properties of the Riemann tensor. 
Thus all of the contributions to the quantity A(n /i n /i ) 2 / (n^n^) 2 come from the V-component 
of n M . 

Fluctuations in n v indicate that an outgoing null geodesic may take either a longer 
or shorter than average affine time in reaching some distance from the black hole. In this 
sense the results tend to agree with the heuristic picture of Ford and Svaiter [ll| • This would 
correspond to an uncertainty in the knowledge of which wavepacket was under consideration. 
Recall that in the derivation of the Hawking effect the ingoing wavepackets were controlled 
by an integer, n, which allowed for successive ingoing wavepackets. A fluctuation in n v 
would mean that the ordering of these successive wavepackets could be disrupted. This 
would not have any observable effect on the outgoing radiation if the ingoing state is the 
vacuum. But it is also possible to have stimulated emission in addition to the thermal flux. 
This would be caused by particles that are initially present during collapse and has been 



studied by Wald [35J. In this case the stimulated emission only occurs at early times while 
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the late time behavior remains thermal. Fluctuations in n could allow for fluctuations in 
the arrival time of particles originating from this stimulated flux. 

Note that effects near the horizon do not necessarily translate to effects observed at J^ + . 
Indeed, we were restricted to considering fluctuations only along a short segment of the 
outgoing geodesic. In order to discuss observations made by a distant observer, one would 
really have to be able to follow the evolution of the separation vector all the way out to the 
observer. In doing so one would see whether the fluctuations integrate to produce a large 
effect or not. The presence of the effective potential barrier prevents us from tracking the long 
term evolution of the fluctuations and so at this point discussing observations of a distant 
observer is merely speculation. Future research may help resolve this issue. In particular, 
particle creation near a moving mirror has been shown to be analogous to the Hawking effect 
for black holes, with an identical thermal spectrum obtained for specific mirror trajectories 
[H, H3, HH]. Further insight may be obtained into the horizon fluctuations of black holes 
by considering fluctuations in the trajectory of the mirror. The benefit of considering a 
moving mirror is that it lacks the complicated potential barrier of a black hole, allowing one 
to integrate out to an observer. 



VIII. SUMMARY AND CONCLUSIONS 

Deviation of outgoing null geodesies in a fluctuating Schwarzschild geometry have been 
considered. The reference geodesic was taken to be that outgoing null geodesic which gener- 
ates the future horizon in the average (Schwarzschild) space-time. Fluctuations of the space- 
time were induced both actively and passively. The active fluctuations were first described 
by a scalar graviton model, consisting of a conformal transformation of the Schwarzschild 
metric where the conformal factor is a quantized scalar field, and then by a graviton model 
that is constructed as a linear quantum tensor field where the mode functions are taken to 
be the classically allowed even parity Schwarzschild perturbations. The passive fluctuations 
arise from stress tensor fluctuations of a free scalar field. In all models the ingoing pertur- 
bation is taken to occupy a multimode squeezed vacuum state. The modes in question are 
wavepackets being sent into the black hole after the collapse process. 

For all models, fractional fluctuations of the quantity n M n M were calculated. The vacuum 
level fluctuations are very small for large black holes, being of order (£ p /M) 2 for active 
fluctuations and of order (£ p /M) 4 for passive ones. However, the fluctuations could be 
boosted by increasing the amount of squeezing. An upper bound on the squeeze parameter 
was imposed by restricting the energy density of the ingoing field to fall within the allowed 
limits of semiclassical theory ignoring backreaction. This upper bound was then used to 
estimate an order of magnitude for the fractional fluctuations, A(n M fi A1 ) 2 / '(n^n^) 2 . It was 
found that the fractional fluctuations could in principle become of order unity. 

Such large fluctuations would seem at first sight to have a dramatic effect upon the 
outgoing modes which carry the created particles. However, the fluctuations come from the 
^-component of not the [/-component. This implies that the primary effect of these 
fluctuations is on the time delay of individual wavepackets, which would only be observable 
in stimulated emission, but not in spontaneous emission arising when the quantum state 
is the in-vacuum. This suggests that the thermal nature of the Hawking radiation is quite 
robust and is not altered by the type of enhanced fluctuations we have studied. 
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APPENDIX A: SCALAR FIELD NORMALIZATION 

The scalar field normalization is found via the Klein-Gordon norm 

(lpjn,tpfn') = i j \ Tp] n V M Tpj'n^J = S U >6 mm f. (Al) 

The surface over which the integral is performed is naturally J^~, past null infinity, and 
the coordinates appropriate in this region are ingoing Eddington-Finkelstein coordinates 
(v, r, 8, ip). In these coordinates, J'" is a three-surface isomorphic to S 2 x R, or a sphere at 
infinity further parameterized by v. The normal to this surface is therefore in the direction 
of d/dr and the preceding integral is 

/•«>/«-* 

,2jn / j„. / „/.* \-ir 



% dZ„ U* jn W Vj-vJ = i J r 2 dn J dv u* n v r 

The field ipj n is a scalar field, for which V M = d^. The metric for ingoing Eddington- 
Finkelstein coordinates is off diagonal, and the (v,r) sector is 





. A <r=L , (A2) 

From this one finds 

/ 2M" 

d r = g rv d v + g rr d r = d v +(l )d r (A3) 



and may proceed to calculate 



oo 



i r dfl dvib nr , 

' S 2 I 



Ipj'n' ■ (A4) 



The calculations are straightforward, but one must bear in mind that while the advanced 
time is implicitly defined in terms of r, when taking the derivative with respect to v or r, 
the other variable is held fixed. Additionally, the r dependence is the same for both ipj n and 
ipj'n' so the r derivative terms cancel. Examining the norm in the asymptotically flat region 
— > r —>■ oo, Eq. flA4j) becomes 



i 



2n £j 



dnY2 m (e,ip)Y L , m ,{e,ip) 



x 



0'+i)£j rti+tyj , F*(r)F ,(r) , 

duo I dj^-^ie ^^n^yn j dv ^ 9ve -^vy (A5 ) 



From the normalization condition of the spherical harmonics, the integral over the sphere 
gives a product of delta functions, Sw6 mm '. The function F(r) was included in the definition 
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of the ingoing null field to reflect the presence of an effective potential barrier due to the 

space-time curvature, but F(*y~) = f where r — > oo at J?~. Lastly, e luJV d v e~ lul ' v = 
—i(co + oj')e l ^ u ~ w > v so the integration with respect to v gives a delta function 7iS(u — uj'). 
This reduces the normalization condition at infinity to 

(lpj n ,1pj' n ') = (e^Su'Smm' / du = 5 U i5mm>- (A6) 

So as given, ipj n is properly normalized. 



APPENDIX B: GRAVITON MODE NORMALIZATION 

Some care is required when fixing the normalization of the ingoing graviton wavepackets. 
We choose to fix the normalization in the asymptotically flat-space region of — > oo 
by setting the vacuum energy of each mode to \uj. One would like to first calculate the 
effective energy density of the wave, integrate over some volume, and set the result to 
\u. The typical approach in calculating the energy density would be to use Eq. 35.70 in 
Ref. [HI], which relates the effective energy density of a gravitational wave to derivatives 
of the perturbation. The problem is that in the Regge- Wheeler gauge, where calculations 
are most easily performed, ^(j)^ is not well-behaved at infinity. In fact, Hg(r), H\(r), 
and H 2 (r) all grow linearly with r as r* — > oo. While Eq. 35.70 of Ref. [39] is gauge 
invariant, it implicitly involves some averaging and assumptions about covariant divergences 
at infinity that do not hold in the Regge- Wheeler gauge. Alternatively, it should be possible 
to transform the perturbations to the radiation gauge, where they are well-behaved, but 
the calculations become more difficult in the radiation gauge. Instead, we choose to take 
a somewhat circuitous route by first calculating the (gauge invariant) Riemann tensor at 
infinity. We then transform the Riemann tensor to Cartesian coordinates and make use of 
a special relation between the Riemann tensor and gravitational waves in the Transverse 
Tracefree (TT) gauge. Calculating the effective stress tensor of a gravity wave in the TT 
gauge is then straightforward using 

T«// = J^/tfTr ( y ab(TT)y ,^TT x*^ «6(TT)v (m) 
± P-v 32tt \*(3H/iI*(j),i/ I ^ V* (j)a&,fJ *0V ' ^ ' 

Here the brackets indicate a spatial average over several wavelengths. The components of 
$m have a simple relationship to components of the Riemann tensor, ^^u. i 00 = —2R koi- 

Since ^(j)^ has a simple sinusoidal time dependence, it follows that ^f^ui o = ~ iLJ ^fj)ki anc ^ 

so ^fj) k i = 2u;~ 2 -Rofco«- In the asymptotically flat space at infinity, 

Let be the transformation matrix from spherical to Cartesian coordinates. Then in 
Cartesian coordinates 

RqP'Ou' = ^0 Ra'P'ii'v' = ^0 Jl'Jp'Jp'Jv'RjXpcr = J pi J°> Ro\Oa (B3) 
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while the TT metric perturbation is \1/^ T 
calculations, one finds 



uj 2 Jf-RoAOo-- Carrying out the above 



(.m 



T e ^ = e- 1 (P/(cos£) {l + (1 - £(£ + l)) 2 ) 



2cos6Pt(cos6) (sin 2 0P/(cos0) - cos#P;(cos#))) 

,CG'uu' 



X du du 



167rr 5 



cos((r + t + 5)(u;-u/)). (B4) 



Next, for ingoing null waves, the total energy on the surface of a sphere is the same as the 
flux through the sphere. The flux through the surface of a sphere is therefore 

$ = f T^ f r 2 dn = e _1 27r / sin0d0 (P/(cos0) (l + (1 - £(£ + l)) 2 ) - 
2cos#P;(cos#) (sin 2 0P/(cos0) -cos0P;(cos0))) 

,ccw 



x I du I du 



16n 



cos({r + t+8)(u-u')) . (B5) 



Consider first the angular integral. Let x = cos 9. The defining differential equation for 
Legendre Polynomials is 



(1 - x 2 )P"{x) - 2xP £ '(x) + £{£ + l)Pi(x) = . 
and the associated Legendre functions may be defined by 



P 



(m) 



(-ir(i-x 2 r/ 2 —- p(x). 



They satisfy the orthonormality conditions 



dxPp(x)P t 



dx r 



2{£ + m)\ 



and 



1 , P £ {n) (x)P, M 
1 1-x 2 



(2£+l)(£-m)\ 
0, m n 



(B6) 



(B7) 



(B8) 



m = 7^ • 



(l+m)! 
m(£—m) \ ' 

oo, m = n = 



(B9) 



Using Eq. ( 1B6I) to rewrite P"(x) in terms of P^(x) and P^(x), the angular integral in Eq. ( 1B5I) 
becomes 



P/(x) (1 + (1 - £(£ + l)) 2 ) - 2x 2 (P;(x)) 2 + 2£{£ + l)xP^x)P e (x) 



(BIO) 



The first term satisfies the equal £ orthonormality condition, Eq. (1B8|) . so 



(! + (!-£(£+ 1))) / dxP 2 {x) = (! + (!- £(£ + 1))) 



2£ + l 



(BU) 



27 



The second term is a little more complicated. Begin by integrating by parts 

• 1 r l 



dxx 2 {P' l {x)) 2 =[x 2 P' l {x)PfXx)\ l _ 1 - / dx[2xPfXx)P' i {x)+x 2 P t {x)P ! ;{x)]. (B12) 

J -i 

Using Eq. (jB6|) . x 2 P"(x) = P"{x) — 2xP' l [x) + £(£ + l)Pi{x) and the preceding is now 

[x^P'^P^W -£(£+ 1) y dxPt(x)P e (x) - j dxP e (x)P"(x). (B13) 

The second term here satisfies the orthonormality condition, while integrating the third term 
by parts and using Eq. (]B7j) gives 



r 2p , npn ,l [p ,, , pf ,1 21(1+1) /- 1 P, (1) Qr)P, (1) (:r) 



(B14) 



where the first and second terms cancel. Meanwhile, the last term satisfies Eq. (1B9|) and 
finally 



2,„„ u2 (« + !) ! 2f(f + 1) 



.^(w^-w- (B15) 



Consider now the last term of fIBlOl) . An integration by parts gives 

dxxPl(x)Pe(x) = [xP e (x)Pi(x)\ l _ 1 - [ dx xP' e {x)P t {x) - [ dxP e (x)P £ (x). (B16) 
i J-i J-i 

The second term on the right is the same as the original integral, while the third integral 
on the right satisfies Eq. (1B8j) . Furthermore, Pe(l) = 1 while P|(— 1) = (— 1) so that 
[xPi{x)Pi{x) ] \_ l = 2 and therefore 

1 2£ 
dx xP^(x)P t (x) = ——. (B17) 



i 



2£ + 1 



Adding up all the results gives 

sin 6d6 (P/(cos 6) (1 + (1 - £(£ + l)) 2 ) 



- 2cos#P;(cosfl) (sin 2 0P/(cos6>) - cos0P;(cos0))) 

= ^[2 + ^+l)(£ + 3)]-|^ (B18) 

The flux through the surface of the sphere is now 

® = £- 1 L j ' du j dJ CG '^' cos ((r + t + 6){u - J)) (B19) 

where 
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The total energy of the wavepacket may now be found by integrating the flux through the 
surface of the sphere for all time, J dt$>. Writing the cosine function in terms of exponentials 
it is clear that f dtcos ((r + t + 5)(oj — lu')) = 2tt5(oj — oj') cos((r + 5){oj — to'))- The delta 
function takes care of the integration over u', and setting the total energy to uj/2 (h = 1) 
gives the condition 

e -i£ / dJ^- = -00. (B21) 



2 2 

For a wavepacket sharply peaked in frequency, we may approximate 



uTdcj « uTAu; = uj 2 e, (B22) 



then 

A, m (u) = C= -±=. (B23) 
V txLuj 

This normalization constant is the one contained in the radial function Z(r), so that 
when the normalized wavepacket is written 



V+^i du 



= / ~i==^K> (B24) 



it is understood that the radial function takes the value -ZYr* — > oo) 



APPENDIX C: SQUEEZED STATES 

A squeezed state is the natural state for a quantum mechanically created particle occu- 
pying an in-vacuum state represented in an out-Fock space. Squeezed quantum states are 
generated via the unitary displacement and squeeze operators. This Appendix provides a 
brief summary of the relevant ideas and results for squeezed states, primarily following the 
notation found in Ref . [13] ; see also Refs. IH, S3] • 

Squeezed states are generate using the unitary displacement and squeeze operators. The 
displacement operator is 

D(a) = e aaWa . (CI) 
One can check that D(pi) transforms a and a) as 

D ] (a)dD(a) = a + a and D ] (a)d ] D(a) = a) + a*. (C2) 

The squeeze operator is 

S(C) = exp [\Ca 2 - |C(a f ) 2 ] , C = pe W - (C3) 

Here the squeezing parameter, (, is an arbitrary complex number. One may show that the 
squeeze operator transforms a and a) as 

S\QdS(Q = acoshp-a^^sinhp (C4a) 

and 

5 f (()tfS{Q = ^coshp - ae- ie sinhp. (C4b) 
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1. Multimode Squeezed State 

Suppose the state |0, () is a multimode squeezed state of the form 

m 

io,o=n^)i°>- ( c5 ) 

The expectation values under consideration require the calculation of terms such as 
J2j Sfc(C) 0\djd k \0, C), J2j Zlfc(C 0|aja fc a r a s |0, C), etc. Consider for example 

m m 

(C,0|a}a fc |0,C> = (0\l[s\(e)a]a k l[S(( e )\0). (C6) 

This expands as 

(o|(5t(c m )5t(c m _i)...5t(c„)a j 5(c„)5(c rH -i)-5(U) 

X (5 t (Cj5 t (Cm-l)...^ t (Cn)a^(Cn)5(Cn + l)..^(Cm)) |0). (C7) 

The action of S(Q) on aj and Oj is 

S^COajS^O) = (a] coshp^ - di sinh p e )5 je + a](l - (C8a) 

and 

^(O^SXOO = («j coshp^ - a] sinh p^t^ + a,(l - o^). (C8b) 
The action of a range of squeezing is then 

S ,t (Cm)-5 ,t (C„)a}5'(Cn)...5 , (Cm) = (aj cosh p, - a,- sinh Pj )G nm ti) + - Q nm {j)) , (C9) 
where the integer step function is defined as 

1, n < j < m, 

QnrnU) = < (CIO) 

0, otherwise. 

It follows that 

^(C,0|a]a fe |0,C) = 5^(C,0| (ajcoshpj - a,- sinhp^e^^') + aj(l - 9 nm (j)) 

jk jk 

x [(a* cosh p k - d{ sinh p k )Q nm (k) + a fc (l - 9„ m (A;)) |0,C) 

= ^2 l(~®nm(j) sinhpj) (-6 nm (fc) sinhp fe )] <5 jfc . (Cll) 



To summarize, one finds 



5^(C,0|ajafc|0,C) = ^2i( 1 + <3> nm(j) (cosh pj - 1)) (-6 nm sinhp fc )] 5 jk , (C12a) 
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0|a]4|0, = J2 K- e n m {j) sinh Pj ) (1 + e nm (k) coshpfe)] (J ifc , (C12b) 

jk jk 

y](C,0loj4l°»0 = X] [(- e nm(j) sinhpj) (-0 nTO sinhp fe )] 5 jfe , (C12c) 

and 

^(C, 0|a}a*|0, C) = J] [(1 + e nm (j)(cosh Pj - 1)) (1 + Q nm {k) coshp fe )] 6 jk . (C12d) 

Note that the result of Eq. flC12dj) contains a 5^ that makes the sum divergent. Renormal- 
ization is therefore taken to correspond to restricting the sum over modes to those occupying 
an excited squeezed state mode, i.e. Y^jL n - For products of two operators this is equivalent 
to normal ordering; equations flC12aj) - flC12dl) become 

yj(C) 0|ajcifc|0, = X] ~~ coshpj smb. p k 5j k , (C13a) 

jk jk 

0\a\al\0, C) = ^ - sinh pj cosh p k S jk , (C13b) 

jk jk 

/]((, 0|%afc|Q, C) = ^ sinh p i sinh p k 5 jk , (C13c) 

jk jk 

and 

X(C,0|ata fe |0,C) = ^ cosh p 3 - cosh pfc^-fc. (C13d) 

jk jk 

The study of scalar field stress tensor induced fluctuations further requires the use of the 
expectation value of four-operator products such as 

Y Y Y XX °\^katal\0, = 

j k r s 

X(°l coshpj - at sinhp 3 )e nm (j) + aj(l - 6 nm (j)) 

jkrs 

a k coshp fc - a\ sinh p fc )6 nm (A;) + a fc (l - Q nm (k)) 
x [(aj coshp r - a r sinhp r )6„ m (r) + a r f (l - nm (r))] 

x [(a* cosh p s - a s sinh p s )6 nm (s) + aj(l - nm (s))] |0) (C14a) 

X K 1 + Onm(j)( COsh Pj - + Q nrn(k) (cosh p k - 1)) 

jkrs 

x (1 + nm (r)(coshp r - 1))(1 + nm (s)(coshp s - l)){S jr S ks + 5 js 5 kr )] 
+ [(1 + 6 nm (j)(coshpj - l))(-sinhp fc nm (A;)) 

x (1 + nm (r)(coshp r -!))(- sinh p s Q n m{s))5 jk 5 rs ] . (C14b) 



x 
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In this case there is the term ((Sj r 8k s + Sj S S kr ) which makes the sum over modes divergent. 
Again, renormalization corresponds to restricting the sums over modes to those modes which 
lie in the range of squeezing. As discussed in in Sect. IVI} restricting the sum over modes 
this way for the four-operator products corresponds, in the limit pi 1, to retaining only 
the fully normal ordered term. The required results are presented here, incorporating the 
restriction on the sum over modes, but omitting the summation symbol for notational sim- 
plification. 

((, 0|djdfcd r d s |0, = cosh pj cosh p k sinh p r sinh p s {Sj r 5 ks + 8j S 6kr) 

+ cosh pj sinh p k cosh p r sinh p s (Sj k 5 rs ), (C15a) 



((, 0\aja k alal\0, () = cosh pj cosh p k cosh p r cosh p s (Sj r 5 ks + Sj S S kr ) 

+ cosh pj sinh p k cosh p r sinh p s (SjkS rs ), (C15b) 



(C, 0\ajdka r al\0, () = — cosh pj cosh p k sinh p r cosh p s (5j r 5 ks + Sj S S kr ) 

— cosh pj sinh p k cosh p r cosh p s (5j k 5 rs ), (C15c) 



(C, 0\djdkdld s \0, — — cosh pj cosh p k cosh p r sinh p s (Sj r 5 ks + Sj S S kr ) 

— cosh pj sinh p k sinh p r sinh p s (5j k 5 rs ), (C15d) 



(C, 0|d^aj.d r a s |0, = sinh p.,- sinh p fc sinh p r sinh p s (Sj r 5 ks + 5j S 5 kr ) 

+ sinh pj cosh p^ cosh p r sinh p s (Sj k S rs ), (C15e) 



((, 0|d^a[,djaj|0, () = sinh p^ sinh p k cosh p r cosh p s {8j r 5 ks + <5 JS (5fc r ) 

+ sinh pj cosh p k sinh p r cosh p s (Sj k 5 rs ), (C15f ) 



(C, 0\djdld r dl\0, C) = - sinh p.,- sinh p k sinh p r cosh p s (5j r 5 ks + 5j S 5 kr ) 

— sinh pj cosh p^ cosh p r cosh p s (Sj k 6 rs ), (C15g) 



(£, 0|a^a{,aja s |0, C) = — s i nn Pj s i nn Pfc cosh p r sinh p s (Sj r 5 ks + Sj S S kr ) 

— sinh pj cosh p fc sinh p r sinh p s (Sj k S rs ), (C15h) 
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((, 0\aja\.a r a s \0, = — cosh pj sinh p k sinh p r sinh p s (Sj r 5 ks + 5j S 5 kr ) 

— cosh pj cosh pfc cosh p r sinh p s {5j k 5 rs ), (C15i) 



((, 0\ajd k alal\0, () = — cosh p.,- sinh p k cosh p r cosh p s (5j r 5 ks + Sj S $, 



kr I 



cosh pj cosh pfc sinh p r cosh p s (Sj k S rs ), (C15j) 



(C, 0\a.jala r al\0, () = cosh pj sinh pk sinh p r cosh p s (Sj r S ks + Sj S 5 



kr I 



+ cosh pj cosh pfe cosh p r cosh p s (Sj k 5 rs ), (C15k) 

((, 0\djdldla s \0, () = cosh pj sinh p fc cosh p r sinh p s (5j r 5 ks + S js 5 kr ) 

+ cosh pj cosh pfc sinh p r sinh p s {Sj k S rs ), (C151) 

(C, 0|a]cifca r a s |0, () = - sinh pj cosh p fc sinh p r sinh p s (5 jr 5 ks + ^s^fcr) 

— sinh pj sinh pfc cosh p r sinh p s (Sj k S rs ), (C15m) 

(C, 0|a]afcajcig|0, () = — sinh p^ cosh p k cosh p r cosh p s (Sj r S ks + ^js^fcr) 

— sinh pj sinh p k sinh p r cosh p s (<5jfc5 rs ) , (C15n) 

(C, 0|ata.fca r aj |0, () = sinh p.,- cosh p k sinh p r cosh p s (Sj r 5 ks + Sj S 5 kr ) 

+ sinh pj sinh p k cosh p r cosh p s (8j k 5 rs ), (C15o) 

and 

((, 0|aja fc aja s |0, () = sinh pj cosh p k cosh p r sinh p s (5j r 5 ks + 5j S 5 kr ) 

+ sinh pj/c sinh p k sinh p r sinh p s {5j k 5 rs ). (C15p) 



[1] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[2] S. B. Giddings, Phys. Rev. D74, 106005 (2006), hep-th/0605196. 

33 



W. G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

W. G. Unruh, Phys. Rev. D51, 2827 (1995). 

T. Jacobson, Phys. Rev. D44, 1731 (1991). 

T. Jacobson, Phys. Rev. D48, 728 (1993), hep-th/9303103. 

T. Jacobson, Phys. Rev. D53, 7082 (1996), hep-th/9601064. 

S. Corley and T. Jacobson, Phys. Rev. D54, 1568 (1996), hep-th/9601073. 

T. Jacobson and D. Mattingly, Phys. Rev. D63, 041502 (2001), hep-th/0009052. 

J. W. York, Phys. Rev. D28, 2929 (1983). 

L. H. Ford and N. F. Svaiter, Phys. Rev. D56, 2226 (1997), gr-qc/9704050. 
R. Parentani, Phys. Rev. D63, 041503 (2001), gr-qc/0009011. 

C. Barrabes, V. P. Frolov, and R. Parentani, Phys. Rev. D62, 044020 (2000), gr-qc/0001102. 

B. L. Hu, A. Raval, and S. Sinha (1999), gr-qc/9901010. 

K. L. Tuchin, Nucl. Phys. B553, 333 (1999), gr-qc/9807065. 

C. H. Wu and L. H. Ford, Phys. Rev. D60, 104013 (1999), gr-qc/9905012. 

B. L. Hu and A. Roura, Phys. Rev. D76, 124018 (2007), arXiv:0708.3046 [gr-qc]. 

S. B. Giddings, Phys. Rev. D76, 064027 (2007), hep-th/0703116. 

R. D. Sorkin and D. Sudarsky, Class. Quant. Grav. 16, 3835 (1999), gr-qc/9902051. 

N. G. Phillips and B. L. Hu, Phys. Rev. D62, 084017 (2000), gr-qc/0005133. 

J. Borgman and L. H. Ford, Phys. Rev. D70, 064032 (2004), gr-qc/0307043. 

L. H. Ford and C.-H. Wu, Int. J. Theor. Phys. 42, 15 (2003), gr-qc/0102063. 

B. L. Hu and E. Verdaguer, Living Rev. Rel. 7, 3 (2004), gr-qc/0307032. 

L. H. Ford and T. A. Roman, Phys. Rev. D72, 105010 (2005), gr-qc/0506026. 

R. T. Thompson and L. H. Ford, Phys. Rev. D74, 024012 (2006), gr-qc/0601137. 

R. M. Wald, General Relativity (University of Chicago Press, 1984). 

T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

L. A. Edelstein and C. V. Vishveshwara, Phys. Rev. D 1, 3514 (1970). 

C. V. Vishveshwara, Phys. Rev. D 1, 2870 (1970). 
F. J. Zerilli, Phys. Rev. Lett. 24, 737 (1970). 

F. J. Zerilli, Phys. Rev. D 2, 2141 (1970). 
R. H. Price, Phys. Rev. D 5, 2419 (1972). 
R. H. Price, Phys. Rev. D 5, 2439 (1972). 



34 



[34] S. A. Teukolsky, Phys. Rev. Lett. 29, 1114 (1972). 
[35] R. M. Wald, Phys. Rev. D13, 3176 (1976). 

[36] P. C. W. Davies and S. A. Fulling, Proc. Roy. Soc. Lond. A348, 393 (1976). 
[37] P. C. W. Davies and S. A. Fulling, Proc. Roy. Soc. Lond. A356, 237 (1977). 
[38] L. H. Ford and A. Vilenkin, Phys. Rev. D25, 2569 (1982). 

[39] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (W.H. Freeman and Co., San 

Francisco, 1973). 
[40] C. M. Caves, Phys. Rev. D23, 1693 (1981). 
[41] D. Stoler, Phys. Rev. Dl, 3217 (1970). 
[42] D. Stoler, Phys. Rev. D4, 1925 (1971). 
[43] H. P. Yuen, Phys. Rev. A13, 2226 (1976). 



35 



